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ABSTRACT 

In this work we present a detailed analysis of the global and fine structure of 
four middle-mass disc galaxies obtained from hydrodynamic simulations in a ACDM 
scenario. These objects have photometric disc-to-total ratios in good agreement with 
those observed for late- type spirals, as well as kinematic properties in agreement with 
the observational I-band Tully-Fisher relation. We identify the different dynamical 
components at zero redshift on the basis of both orbital parameters and the binding 
energy of stars in the galaxy. In this way, we recognize a slowly rotating centrally con- 
centrated spheroid, and two disc components supported by rotation: a thin disc with 
stars in nearly circular orbits, and a thick disc with orbital parameters transitional 
between the thin disc and the spheroid. The spheroidal component is composed mainly 
by old, metal-poor and a-enhanced stars. The distribution of metals in this component 
shows, however, a clear bimodality with a low-metallicity peak, which could be related 
to a classical bulge formed from rapid collapse at early times, and a high-metallicity 
peak, which could be related to a pseudo-bulge formed from instabilities of the inner 
disc. The thin disc appears in our simulations as the youngest and most metal-rich 
component, with median stellar ages ranging from 3.8 to 6.7 Gyr. The radial distri- 
bution of ages and colours in this component are U-shaped: the new stars are forming 
in the inner regions, where the galaxy is bluer, and then migrate through secular pro- 
cesses reaching the outer parts. Finally, we also find, in all simulated galaxies, a thick 
disc containing about 16 per cent of the total stellar mass and with properties that 
are intermediate between those of the thin disc and the spheroid. Its low-metallicity 
stars are a-enhanced when compared to thin disc stars of the same metallicity. The 
structural parameters (e.g., the scale height) of the simulated thick discs suggest that 
such a component could result from the combination of different thickening mecha- 
nisms that include merger-driven processes, but also long-lived internal perturbations 
of the thin disc. 
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1 INTRODUCTION 

Spiral galaxies have long been recognized to contain two 
main stellar components: a dynamically hot spheroidal com- 
ponent composed by old and metal-poor stars, and a cold 
disc of young and more metal-rich stars. 

In the past three decades, a wide range of observa- 
tions has revealed the existence in almost all resolved spi- 
ral galaxies of yet another component, a thick stellar disc. 
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This component wa s originally detected in the Milky Way 
using star counts (jGilmore ReidI [1983 ) and as an ex- 
cess r ed flux at large galactic latitud es in exte rnal galaxies 
(Dalc anton B ernstein 2002; Yoach im &; Dalc anton 20061). 
It is now k nown to be an ubiquitous component of sp i- 
ral galaxies toalca nton k Bernsteinll2002l : ISeth et al.ll2005l l. 
but for many years it remained unclear whether it was a 
truly distinct component of the Milky Way thin disc. Evi- 
dences from their distinct nature came from detailed kine- 
matic and chemical abundance studies (e.g. Fuhrma nnll 19981 : 
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Soubiran et alll2003l : iBensbv et al1l2003l . l2005l : iReddv et all 
2QQ6I ). 

In general thick discs are characterized by an 
exponential- fit scale height, Zh, larger than the th i n disc 
by a factor of ^2 (e.g. IPalcanton BernsteinI l2002l : 
lYoachim fc Dalcanton 2006). In our Galaxy and in NGC 
891, a Milky Way analogue, the Zh value seems to be 
roughly constant with radius { Robin et al.ll2003l : llbata et"al] 
I2OO9I I. Relative to the galactic thin disc, thick disc stars 
have both larger veloc ity dispersions and slower net rotation 
(|Soubiran et al.ll2003l ). It conta ins 10 to 35 per cent of the to- 
tal baryonic mass of a galaxy (jYoachim Dalcanton|[2006l : 
ICollins et aDl2010l \ and appears to be anti correlated with 
the bulge mass: galaxies with small bulges have larger thick 
discs. Concerning their stellar population, thick disc stars 
are older and more meta l -poor than stars in the thin disc 
(e.g. iGilmore et aD Il995l : \Fuh rman ^ Il998l ). They are also 
significantly enhanced in a-elements com pared to thin disc 
stars of com parable iron abundance (e.g. Fuhrmann "l998|; 
iBensbv et al. 2003, 2005; Reddv et al. 2006), suggestinga 
star formation tim e-scale shorter than 1 Gyr (iFuhrmannI 
Il998l : iRuchti et alllioiQ ) . 

Several formation mechanisms have been proposed in 
the literature to explain the origin of the different galac- 
tic components. In the prevailing picture of galaxy forma- 
tion, haloes of dark matter acquire angular momentum via 
tidal torques from interacting structures. As gas cools and 
condenses into their central parts, star-formin g galaxies ap- 
pear with a thin disc structure ([Fall Efstathioul I198QI : 
IWhitelll984l ). Hydrodynamical simulations, however, have 
shown that cold flows play a fundamenta l role in the build- 
ing up of di scs in low-mass haloes (e.g. 'Keres et al.l l2005l : 
IPekel Bir nboim 2006). The cold gas flows rapidly to the 
centre of galaxies from fllamentary structures around haloes, 
where it may be swiftly accreted into the central galaxy. The 
lower initial gas temperature leads to shorter cooling times 
and to the growth of stellar discs at higher redshifts than 
predicted by the standard model (Brooks et al. 2009). 

Within these pictures, two main formation scenarios 
have been proposed for the bulges of spiral galaxies (see 
iKormendv KennicuttI |2004| . for a review). In classical 
bulges, most stars are formed during the early phase of inten- 
sive star formation following collapse of the proto-galaxy and 
subsequent mergers. In another scenario, boxy or pseudo- 
bulges can also be developed through dynamical instability 
of the inner disc. 

Thin disc instabilities are also considered as the main 
responsible for the thick disc formation. In most of the pro- 
posed scenarios the thick disc stars are formed in the thin 
disc and are dynamically hea ted throu gh either the accre- 
t ion of galaxy satellites (Qui nn et al J' 1993; Walke r et al.l 
[l996l : IVillalobos Helmii 2008) or internal perturbations of 
the thin disc from long-lived mechanisms (e.g., the stel- 
lar scattering via resonant interactions with spiral arms 
iRoskar et al ]|2008': 'Schon rich k BinnevI [2009al lbl) or short- 
lived events (e.g., perturbations d ue to the thin disc clu mpy 
stellar structure at high redshifts iBournaud eralll2009l ). In 
other category of models, the thick disc stars are assumed 
to be formed eith e r in situ during/after gas-rich mergers 
iBrook et al ] 12004 12005'), or in external satellites that are 
disrup ted and incorporated into the thick disc (^Abadi et al.l 
l2003bl ). 



Numerical hydrodynamic simulations have been recog- 
nized as a powerful tool to study galaxy formation and 
evolution in its full cosmological context. The main dif- 
flculty in reproducing realistic disc-dominated galaxies is 
the transfer of gas angular momentum to the dark mat- 
ter (Navarro Whit3 ll994V which results in small simu- 
lated discs with too high bulge-to-disc ratios. This catas- 
trophic loss of angular momentum can be prevented by 
modelling processes (e.g., the feedback from supernova ex- 
plosions) that regulate star formation and the rapid con- 
sumption of gas at early times. Although the observed 
large fraction of late type gala xies is still far from be- 
ing reproduced in simulations (|Weinzirl et al.1 l2009l ). m 
the past decade some simulations have been able to re- 
produce ensemble properties of disc galaxie s with gaseous 
discs at lower re solutions ([Dommguez-Tenreiro et al.lll998l : 
ISaiz et al.|[200ll l. as well as to produce individual exam- 
ples of realistic disc galaxies in ACDM universes (e.g. 
' II — " — 1 1 • II 1 1 ' ' ^ I 

Sommer-Larsen et al. 2003; Abadi et al. 2003b; Brook e t al.l 
2004; Governato et al. 2007, 2009; Scannapieco et al. 20Cii). 
These simulations have generated extended stellar disc com- 
ponents with high speciflc angula r momentum and with a 
signiflcant fraction of young stars (jScannapieco et aLll2008l l . 
Although an old and significan t sphe r oid is almost always 
present tSommer-Larsen et al.l [20031: lAbadi et al.l l2003al : 
Ibovernato et al.ll2007l : IScannapieco et al.ll2008l l. these discs 
show realistic exponential scale lengths and integrated 
colours. Moreover, they fall on the region defined by the 
observational I-band Tully-Fisher relation and have bulge- 
to-disc ratios or stellar mass-to-light ratios similar to obser- 
vational counterparts. The presence of a thick disc compo- 
nent is evidenced in some simulated discs by both a double 
exponential vertical structure, and a velocity dispersion ver- 
sus age re lation th at resemble those observed in the Milky 
Way ( Abadi et al.ir2003b : Brook et al. 200 4). 

The aim of this and forthcoming papers is to study the 
formation mechanisms of the different galactic components 
(bulge, halo, thin disc and thick disc) by means of cosmo- 
hydrodynamic simulations. For this purpose we follow the 
evolution of four different disc galaxies in a ACDM cosmol- 
ogy, using an Ope nMP parallel version of the DEVA code 
(iSerna et al ] [2003') which includes the chemical feedback 
and c ooling methods described in iMartmez- Serrano et al.l 
(2008). The present work is focussed on their global and fine 
structure at z = and their comparison with observations. 
We will show in this paper that the dynamical and chemical 
properties of these simulated galaxies have values consistent 
with observations and contain valuable information about 
their formation processes. 

This paper is organized as follows. In Section 2 we 
briefiy describe the simulation code and the main character- 
istics of our simulations. In Section 3 we present the general 
properties of our simulated galaxies at 2; = and compare 
with observational disc properties. In Section 4 we identify 
the dynamical components in the simulated galaxies (sub- 
section 4.1), to perform a detailed description and compar- 
ison of these components with observational references. We 
analyse the vertical density profiles (Section 5), stellar ages 
(subsection 6.1), age and colour gradients (subsection 6.2) 
and chemical properties (Section 7) of such components. Fi- 
nally, in Section 8 we summarise and discuss our main re- 
sults. 
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2 THE SIMULATIONS 

All simulations were performed with an OpenMP parallel 
version of DEVA, an AP3M+SPH code specially designed 
so that conservation laws (e.g. momentum, energy, angu- 
lar momentum and en tropy) hold as accurately as possible 
(see ISerna et al.l l2003l . for details) . The code also includes 
methods to account for subresolution physical processes such 
as star formation, gas restitution from stars, chemical feed- 
b ack and cooling. These method s are extensively described 
in iMartmez- Serrano et al. We will summarize them 

below. 



2.1 Star formation 

Star formation (SF) is implemented through 
Schmidt- like law: 



dt 



Kennicutt- 



(1) 



where p is the gas density, c* is a dimensionless SF efficiency 
parameter, and tg is a characteristic time-scale chosen to 
be equal to the maximum of the local dynamical time and 
the local cooling time. We consider that gas particles are 
eligible to form stars if they are located in a region with 
a convergent flow and a gas density higher than a density 
threshold . At each time-step we compute the probability 
that those eligible gas particles form stars in the considered 
time interval, and draw random numbers to decide which 
of them will form stars. Then, each of these randomly se- 
lected gas particles is transformed into one stellar particle, 
which is assumed to be a single stellar population (SSP) of 
a given age and mass (see Table [1]). In the simulations pre- 
sented in this paper, stellar masses within each SSP were 
assumed to be distrib uted according to the Salpeter initial 
mass function (IMF) (|Salpeterlll955l ), with a mass range of 
[M^,M,] = [O.1,1OO]M0. 

Since our cosmological simulations do not resolve star- 
forming molecular clouds, in the above scheme the choice 
of the yo* value is based on the fact that the ISM can be 
represented by a lognorm al density probability distribution 
function (PDF) (see, e.g., Kravtsov 2003: Wada N ormanI 
I2OO7I I. The amount of gas eligible for star formation is rep- 
resented by the high density part of the PDF, which in turn 
is a function of the total galaxy gas mass and turbulence. 
The choice of the value must allow for the high-density 
star- forming part of the PDF to be well resolved or, at least, 
to contain a converged amount of star- forming ma ss for the 
adopted numerical resolution (jAgertz et al.ll20l"ll ). Despite 
this theoretical dependence of the SF parameters on the lo- 
cal properties, and probably also on redshift, for a given 
simulation such parameters are usually taken as constants. 

Discrete energy injection processes have been found to 
play an important role in the self-regulation of the SF pro- 
cesses in galaxies. Positive feedback effects can help to halt 
the overproduction of stellar mass and, when implemented 
in simulations, help to properly reproduce some proper- 
ties of disc galaxies. The most widely employed model in 
cosmological hydrodynamical simulations assumes a strong 
exp licit feedback combined with a high SF efficiency (e.g. 
IScannapieco et al. 2008 . 2011a: but see also the discussion 
m iKeres et al.1l2009l ). However, a number of authors (see 



'Agert z et al.|[201ll . for a recent work) have shown that the 
sub-parsec processes of SF self-regulation in simulated disc 
galaxies can also be implicitly modelled by means of a low 
SF efficiency. Since our aim in this paper is to test the min- 
imal conditions for the formation of realistic discs, we have 
adopted the latter model for subresolution physics. Accord- 
ing with this approach, we have chosen inefficient SF param- 
eters (see Table [U, which implicitly account for the regula- 
tion of star formation by discrete energy injection processes. 



2.2 Chemical feedback and Cooling 

Methods for chemical feedback and cool ing are also those 
described in lMartmez-Serrano et aL These methods 

account for the full dependence of radiative cooling on the 
detailed composition of the gas, through a fast algorithm 
based on a metallicity parameter, C{T), which takes into 
account the weight of the different elements on the total 
cooling function. Such a method for cooling is important 
to prevent an overestimation of the metallicity-dependent 
cooling rate. 

The release to the ISM of newly produced elements, as 
by-products of stellar evolution and death, is included in 
our code through a stochastic algorithm that completely re- 
moves the assumption of instantaneous recycling of stellar 
ejections. In this algorithm, the metal ejection from stars 
is progressively released to the ISM as increments in the 
metal content of nearby gas particles. The delayed gas resti- 
tution from stars is also taken into account by means of a 
probabilistic approach that reduces the statistical noise and, 
hence, that allows for the study of the inner chemical struc- 
ture of objects with moderately high numbers of particles. 
Furthermore, once the stellar production of metals has been 
released to the ISM, the chemical elements are distributed 
and mixed through the gas component by means of a diffu- 
sion law that models, without the explicit inclusion of stellar 
winds, the subresolution turbulent mixing of metals. 

In the computation of the metal production, our code 
accounts for the full dependence on the detailed chemical 
compositi on of stellar particles by means of the Qij formal- 
ism ( Talbot Arnettl ll973). which relates each nucleosyn- 
thetic product to its sources. Such a Qij formalism is impor- 
tant to prevent a significant underestimation of the [a/Fe] 
ratio in simulated galaxies. In this paper, we hav e used the 
stellar evolution data described in lGavilan et al.l (l2005ll for 
low a nd intermediate mass stars, and IWooslev &: Weaveil 
(|l995l l for high -mass stars. SNIa rates were computed ac- 
cording to iRuiz-Lapuente et al. ' ('2000') and their element 
production according to [iwamoto et al. (1999). We have 
considered the evolution of the following elements: H, ^He, 
i^C, i^C, i*N, i^o, 20Ne, 24Mg, ^^Si, ^^S, ^Ca and ^^Fe. 



2.3 Simulated disc galaxies 

The simulation of each disc galaxy consists of a cosmologi- 
cal zoom-in that includes high-resolution gas and dark mat- 
ter for the flow converging region that generates the main 
object. The rest of the simulation box is sampled by low- 
resolution dark matter particles that account for tidal forces 
over the flow converging region. 

As a flrst step, we performed two low-resolution sim- 
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Table 1. Principal characteristics of simulations: mass of a single baryonic particle (m^ar)? star-formation parameters, density threshold 
(p*) and efficiency (c*), virial radius (rvir), spin parameter of the dark matter halo (Adm), the total (Mvir), stellar (M^^^^^) and gas 
(Mg^g) masses inside the virial radius. 



Galaxy 


mbar 


p* 




C* 


rvir 


Adm 


Mvir 


M^iL 


M^i^s 




[Mo] 


[g cm" 






[/i-i kpc] 




[lOii Mo] 


[lOlO Mo] 


[lOio Mo] 


LD-5003A 


3.82x10^ 


6x10" 


-26 


0.010 


106.12 


0.019720 


1.95 


1.79 


0.63 


LD-5101A 


3.79x10^ 


6x10" 


-26 


0.010 


90.55 


0.034233 


1.21 


1.49 


0.48 


HD-5004A 


3.94x10^ 


6x10- 


-26 


0.010 


127.00 


0.023273 


3.35 


3.74 


1.02 


HD-5103A 


3.78x10^ 


6x10- 


-26 


0.010 


141.62 


0.098032 


4.66 


3.42 


0.88 


HD-5004B 


3.94x10^ 


1x10- 


-25 


0.008 


128.10 


0.023387 


3.44 


3.57 


1.31 


HD-5103B 


3.78x10^ 


1x10- 


-25 


0.008 


140.53 


0.096513 


4.55 


3.36 


0.89 


HD-5004L 


3.12x10^ 


3x10- 


-26 


0.010 


128.47 


0.027249 


3.48 


3.87 


0.96 



ulations that consisted of two different Montecarlo realisa- 
tions of the same cosmological model: a flat ACDM universe 
(with Qa 0.723, Qm 0.277, Qb 0.04, as 0.811 and 
h = 0.7) within a periodic box of 10 Mpc per 

From each of these low-resolution simulations, we se- 
lected two gas-rich objects with a prominent gas disc at 
z = not apparently perturbed from recent major merg- 
ers: the most massive one from a sparsely populated re- 
gion and the most massive one from a crowded region. 
These selection criteria then focus on massive regular spiral 
galaxies extracted from both low and high density regions 
and, hence, with different histories of mergers and accretion 
events. For each selected object, we traced back the parti- 
cles inside its virial radius until the initial redshift 2;init = 31. 



^ Although this box side implies a lack of very massive objects 
and galaxy groups, it has been prove d that it has little ef fect on 
the internal properties of the haloes dPower Knebdl200& ) and 
on the formation of individual galaxies. 



We then computed the convex hull ([Barber et al.lll99"6h en- 
closing these particles at 2;init and substituted all the par- 
ticles inside the convex hull by their high resolution coun- 
terparts. Gas particles outside the hull were eliminated and 
their masses added to the low-resolution dark matter com- 
ponent, thus obtaining the initial conditions for each high- 
resolution simulation. The galaxies obtained in this way have 
a full history of mergers and accretion events in a cosmolog- 
ical context, without any assumptions made for their initial 
conditions beyo nd the c osmology and the initial conditions 
generator used (IPrunet et al. 2008). Table [1] summarises the 
main characteristics of each high-resolution simulation, with 
a gravitational softening of Cg — 400 pc and a minimum 
hydrodynamical smoothing length half this value. In this ta- 
ble, LD-5003A and LD-5101A correspond to objects selected 
in low density (LD) regions, while HD-5004A and HD-5103A 
correspond to objects from high density (HD) regions. All 
these objects (labelled as A) were simulated with identi- 
cal SF parameters. The two objects selected from crowded 
regions experienced major mergers at moderately low red- 
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shifts. In principle, the abihty of their discs to regenerate 
after their last major merger, as well as the resulting prop- 
erties of their fine structure, could be sensitive to the choice 
for SF efficiency. In order to test the possible effects of SF 
efficiency, these two objects were re-simulated with less ef- 
ficient SF parameters (simulations labelled as B: HD-5004B 
and HD-5103B). 

Table [1] also shows the total, stellar and gas masses en- 
closed by the virial radius of each object. We can see from 
this table that our simulated galaxies have virial masses in 
the range 1.2 to 4.7x10^^ M© and stellar masses of 10^°-^, 
placing them among middle-mass spiral galaxies. There- 
fore, they are less massive than the Milky Way (with virial 
mass of -1.2x10^^ Mo) and their masses are comparable 
to the least massive si mulated galaxies presented in recent 
works (e.g. Governato et al.l 120071 : iScannapieco et al.ll200^ : 
iGuedes et aDl201lh . 



3 GENERAL PROPERTIES OF DISCS 

3.1 Bulge-disc decomposition: D/T ratios 

Fig. [1] shows the face-on and edge-on synthetic images of 
the four main simulated objects at 2; = 0, where SSP ages, 
masses and metallicities have been converted into luminosi- 
ties using the Bruzual & Chariot ( 2003) modelfl. In all these 
images, a conspicuous disc component can be appreciated 
even for the object HD-5103A, where the disc structure has 
survived a recent major merger at z ^ 0.3. Although the 
synthetic image of HD-5103A does not show a clear evidence 
of such a merger, we will see in Section [H that this event has 
still a strong imprint in several physical properties. 

In order to quantify the importance of each component. 
Fig. [21 shows the face-on r-band luminosity profiles as well 
as the fits obtained from their bulge-disc decomposition (see 
Table[2]for the fitting parameters in different bands). In each 
band we have performed a one dimensional fit which com- 
bines a Sersic law for the inner region and a simple exponen- 
tial profile in the outer parts. This fit provides a very good 
description of the stellar luminosity profile up to a limit of 
~24 mag arcsec"^. In the outer reg ions, where a second ex- 
pone ntial should be employed (see iMartmez-Serrano et al.l 
I2OO9I ). some deviation is instead observed. Such a deviation 
has however little effect on the global fit quality measured 
with the mean square deviation, x^- 

Recent works have used 2D fitti ng packages d eveloped 
for use in observational studies fe.g. lGadottill2008l ) to per- 
form the bulge-disc decomposition. They have shown that 
different results can be obtained when different fitting tech- 
niques are employed, and also that the effects of neglecting 
bars when they are present can alter the resulting bulge- 
to-total and disc-to-total ratios. In order to test the effects 
of such different techniques, we have also perfo rmed fits t o 
face-on synthetic images using GALFIT (Peng et al .11 20 10 ). 
The resulting fits show similar values for the scale length and 
effective radius, somewhat higher values for the Sersic shape 
parameter n, and minor-to-major ratios very close to one. 



^ Dust effects have not been considered in this work, so, in order 
to minimize these effects, photometric studies have been done 
with galaxies seen face-on. 



The resulting D/T ratios were roughly the same as those 
obtained from ID fits, but with a higher uncertainty due to 
the additional parameters (point-spread function, dust pro- 
cessing,...) required to generate the synthetic images. Due to 
the similitude of results and the fact that our objects have a 
high azimuthal symmetry, we have chosen to show here only 
the one-dimensional bulge-disc decompositions. 

The D/T ratios, measured in our objects from their ID 
fits in different photometric bands, are shown in Tabled In 
general, they are in good ag reement with those expected for 
spiral galaxies fe.g. iLaurikainen et al.|[2007l : lOohama et al.l 
2009 ) and are c omparable with those recently obtained i n 
simulations fe.g. lGovernato et al.ll2009l : lGuedes et al ] |201lh . 
Further, the values obtained from Sersic fits to the total 
light distribution, are also in good correspondence with both 
the luminosity-s ize and stel l ar ma ss-size relations for disc 
galaxies given bv lShen et al.l (|2003l ). In o rder to compare our 
results with the compilation reported bv I Graham WorlevI 
(|2008l ). we have also computed the bulge-disc decomposition 
in the K-band. The median values obtained for the shape 
parameter (n=1.46), the effective radius (re = 0.37 kpc) 
and the scale length (rs =3.71 kpc) of our discs are in good 
agreement with those reported for spiral galaxies. 

Although these results must be taken with cau- 
tion^ due to uncertainties introduced by different fitting 
dGado tti 2008) and deco mposition fe.g. lAbadi et al.ll2003bl : 
IScannapieco et al.l l20ld see also the Section 14.11 below) 
methods, the relatively small size of bulges obtained in these 
simulations is an important result by itself, given the usual 
problems in reprod ucing realistic disc galaxies with co smo- 
logical simulations (jAbadi et al.ll2003al : [BurkertI l2009l ). As 
already quoted in Section [U the main problem hindering the 
formation of large galactic discs is the loss of gas angular mo- 
mentum at early stages of galaxy formation, which results 
in objects with too high bulge-to-disc ratios. This catas- 
trophic loss of gas angular momentum can be prevented by 
including feedback processes avoiding the early in-fall of gas 
and its rapid con sumption (e.g. Sommer-Larsen et al. 2003' 
Governato et al.ll2007l : IScannapieco et al.ll2008l : lAgertz et al] 



201 ih . The cosmological simulations reported in this work, 
where energy feedback is implicitly included in the star for- 
mation parameters, lead to disc galaxies with a bulge com- 
ponent comparable^ to those observed in Sab-Sb galaxies 
(jOohama et al.ll2009l ). 

3.2 Rotation Curves and the Tully-Fisher relation 

Fig. [3] shows the circular velocity curves for the dark mat- 
ter and the baryonic components (dashed and dotted lines, 
respectively) of the simulated galaxies, computed as 



, ^ GM.{< ry 

(r2+ £2)3/2 



(2) 



with i=baryons (bar) and dark matter (dm). Solid lines are 
the circular velocities, V^irci'^) — ^dm(^) + ^barl^)? where r 
stands for the radial distance to the galactic centre. As it can 
be clearly seen from this figure, the baryonic component of 
the galaxy dominates the dynamical properties in the inner- 
most region, while the contribution of dark matter is approx- 
imately constant across the galaxy, dominating the curve in 
the external region. The red solid line stands for the ro- 
tational velocity measured from the gas particle movement. 
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Table 2. Face-on photometric decomposition of simulated galaxies in different bands, re stands for the effective radius of the bulge, r^ 
for the scale length of the disc. Lt and L^, represent the luminosity in the central region and the luminosity of the disc. We also show 
the value of the fit. D/T^ is defined as the disc-to-total ratio obtained using the photometric decomposition, and M as the absolute 
magnitude in each band. 



Galaxy 


Band 


n 












D/T^ 
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We can see that gas particles closely follow the circular veloc- 
ity curve of the objects, except in the innermost region where 
they are too peaky. T his is a common problem in cosmo- 
lo^ical simulations (e.^.lAbadi et aD2003a jGovernato et al.l 



I2OO7I: IStinson et al.H2010l: lAgertz et al.ll2011l : and less severe 
in IScannapieco et al. l2011bl v in which differences between 
the peak circular velocity and the vel ocity at lar^e radii ar e 
larger than in most observed galaxies (Catin ella et al.|[2006[ ) . 

In Fig. |4] we show the I-band Tully-Fisher relation for 
our objects (coloured triangles) as well as for th e sample o f 
spiral galaxies from the catalogue of Spring ob et alj (|2007[ ). 
All the simulated objects lie within the region defined by the 
observational data scatter, although with some systematic 
trend to be fainter than the average for spirals of similar 
rotation speed. 



3.3 Angular momentum 

In order to compare the angular momentum of the different 
components of our simulated galaxies with observations, in 
Fig. [5] we plot the specific angular momentum (j = J/M) 
at 2; = as a function of mass (M) (left panel), and as a 
function of the rotation speed (right pa nel). Recent works 
have made use of similar r epresentations (|Torrev et al.ll201ll : 
IScannapieco et aDl2011al l with the same purpose, since it 
provides a straightforward comparison with observations 
and can help to pinpoint any angular momentum problems. 
In the first plot, the specific angular momentum and mass 
of dark haloes (open squares), stellar bulges (starred sym- 



bols) and baryonic discs (filled squares) are represented sep- 
arately. The mass and angular momentum of the dark mat- 
ter haloes have been estimated using particles within the 
virial radius, while those of the bulges and discs have been 
estimated with star^ within r < Ve and baryons within 
r < Topt = 3.2rs (i.e., in a purely exponential disc ropt 
is the radi us enclosing 83 per cent of the disc mass), re - 
spectively (|Serna et al. 12003: Scannapieco & Tissera 2003). 
From this plot we can see that the disc components and 
the dark haloes have comparable specific angular momen- 
tum, so that disc particles have collapsed conserving, on 
average, their angular momentum. On the contrary stellar 
bulges have been formed from material that had lost most 
of its angular momentum, or that never acquired it. In the 
figure we also include the boxes enclosing the observational 
regions covered by spirals (solid line) and ellipticals (dashed 
line) as given by Fall (1983). 

Given that the observational estimates of specific an- 
gular momentum and mass are derived from observational 
parameters, we also adopt an observational-like approach, 
and compute the specific angular momentum and masses of 
our discs, applying the method of Fall (1983) directly to our 
simulated discs (open circles in left panel of Fig. [5]). The 
specific angular momentum is obtained from the rotation 
curve as j = 2vcrs where Vc = Vcirc(i?25) is the circular ve- 



^ We will refer to stellar particles in the simulation by the generic 
name 'stars', even though they are actually assumed to be single 
stellar populations of a given mass (see Table [T}. 
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Figure 2. Face-on luminosity profiles in the r band (filled tri- 
angles) and their fits obtained with the bulge-disc decomposition 
(red solid lines). 



locity at r25, and masses are derived assuming a constant 
mass- luminosity relation (M/Lb) — 3.0. We see from this 
figure that all simulated discs lie inside the region defined by 
observed spirals and close to the values computed using ob- 
servational methods, although the latter show greater values 
for the specific angular momentum. Similarly stellar bulges 
lie in the region defined by ellipticals, except for some cases 
that still evidence the orbital angular momentum of merging 
structures. 

Another way to quantify any potential angular momen- 
tum problem is to represent j as a fu nction of ^;ro t. This 



app roach has been follo wed recently by Croft et aU (2009l l 
and lAbadi et al.1 (|2003al ) among others. iNavarrol (|l998F pro- 
vides a compilation of specific angular momentum of ob- 
served discs derived using the estimator described above: 
j — 2vcrs. In right panel of Fig. [5] we plot the observed data 
(small dots), along with our objects using three diflFerent es- 
timators for j and Vrot- Filled circles represent the sum of 
specific angular momentum of baryonic particles, and their 
rotational velocity Vrot is that measured at r = 2.2rs. Open 
circles are computed using the observational estimator for 
J, as in right panel of Fig. O and the same estimator for 
'^^rot- Open starred symbols represent the specific angular 
momentum of star particles. Their Vrot is derived using the 
rotation curve of stars only. While all points lie within the 
observationally-defined region, it can be appreciated that 
points derived using the actual baryonic angular momen- 
tum (filled circles) have lower angular momentum than the 
observational- like estimates (open circles). This behaviour 
was also observed in left panel of Fig. \E\ and points to the 
fact that observational estimates, assuming a flat rotation 
curve, can under-estimate or over-estimate the real value of 
the speciflc angular momentum, depending on the trend that 
the rotation curve actually has. The rotation curves of our 
simulated galaxies, decrease with increasing radius owing to 
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Figure 3. Circular velocity profiles of the simulated galaxies. 
Black solid line stands for the total mass distribution (V^irc)? as 
well as green dashed and blue dotted lines for the dark matter 
(Vdm) and baryonic (Vbar) contributions, respectively. The rota- 
tional mean velocities for the gas particles are displayed with a 
red solid line. 
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the relative importance of the bulge, as do most simulated 
discs in the literature. Another effect is that stellar com- 
ponents generally have lower angular momentum than the 
total baryonic component. Given that we take into account 
all the star particles in this estimate, including an old bulge 
and thick disc, this is expected, since these components are 
not fully supported by rotation. 

3.4 Resolution effects 

In order to test for resolution effects in our simulations, we 
have run another instance (HD-5004L) of the HD-5004A ob- 
ject, but with half the spatial resolution and a mass resolu- 
tion of about 8 times lower than in the high-resolution run. 
The gravitational softening length was scaled to Cg = 800 

pc, as well as the density threshold for SF, which was 
set to = 3xl0~^^ g cm~^ in order to roughly have the 
same overall SF efficiency as in HD-5004A (i.e., an object 
stellar to baryonic mass ratio of ~ 85 per cent). The main 
characteristics of this low-resolution simulation, as well as 
some general properties of the resulting object, are shown 
in the last row of Tables [H [2] and [H 

We can see from Table [1] that the resulting virial radius 
and virial mass at ^ = are nearly the same as those of the 
high-resolution counterpart. Both runs give also very close 
values for the stellar (a difference of ^ 4 per cent) and gas 
(a difference of ^ 7 per cent) mass inside the virial radius. 

The photometric decomposition (see Table [2|) also evi- 
dences similar structural properties in the low and high reso- 
lution simulations. The bulge-disc parameters of HD-5004L 
are almost identical to those of HD-5004A (differences al- 
ways smaller than ^ 6 per cent). The only exception is the 
Sersic parameter n, which implies that the stellar distribu- 
tion in the innermost region is less concentrated in the low- 
resolution run. The D/T ratios (see Tabled obtained both 
from photometric fits and a 'kk-means' decomposition are 
also rather similar to those measured in HD-5004A. 

Kinematic properties do not appear significantly af- 
fected by resolution effects either. Indeed, the specific angu- 
lar momenta for the low- resolution dark matter halo, stellar 
bulge and baryonic disc (red symbols in left panel of Fig. 
O lie very close to their high resolution counterparts. In- 
deed, the specific angular momenta for the low-resolution 
baryonic disc (red symbols in right panel of Fig. [5|) lie 
very close to their high resolution counterparts. The cir- 
cular velocity curves of HD-5004L and HD-5004A are also 
quite similar, specially at large radii. As in other works (e.g. 
IScannapieco et aI]|2011bl V the major differences are located 
in the innermost regions, up to 6 kpc, and they are mainly 
due to the different distribution of mass. In particular, the 
peak circular velocity diminishes about 8 per cent in run 
HD-5004L. In this object the contribution of gas to the to- 
tal curve is negligible, so the difference between maximum 
values are mainly due to the stellar distribution. Finally, 
in Fig. [4] we have also included the Tully-Fisher relation 
for the low-resolution counterpart of object HD-5004A, blue 
square dot. This object has a similar total magnitude but 
its baryons are rotating slower, so the final low- resolution 
galaxy is brighter for spirals of similar rotation speed. 

In view of these results, we conclude that the low and 
high resolution objects have converged in their structural 
and kinematic properties for the purpose of our analysis. 
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Figure 6. Mass fraction of stars (black) and gas (magenta) for 
simulated galaxies at z = 0, as a function of ej = jz/Jc- Blue 
dotted, green short dashed and red long dashed lines stand for 
mass fractions of stars in the thin disc, thick disc and spheroid as 
decomposed by 'kk-means' (see text for details). 



4 THE FINE STRUCTURE 

4.1 Dynamical decomposition method 

In order to characterise the properties of different compo- 
nents in the simulated galaxies, we have devised a proce- 
dure to segregate star particles in to the main compon ents, 
follo wing the mode ls proposed bv lAbadi et al.l (|2003bt ) and 
Scannapieco et aP |2009). 

The existence of a thin gaseous disc defines a symmetry 
axis, the 2;-axis, determined by the angular momentum of 
gas particles in the inner regions of haloes at z = 0. We use 
this axis to orientate the object and then compute the eccen- 
tricity parameter, ej, as the ratio of the 2;-component of the 
angular momentum of each baryonic particle, jz, and the an- 
gular momentum expected for a particle with the same bind- 
ing energy in a circular orbit, Jc{E) (Abadi et al. 2003b). In 
Fig.[6]we plot the distribution of this parameter, ej, for stars 
(black line) and gas (magenta line) in all simulated objects. 
As we would expect most gas particles are distributed in 
circular orbits, with ejgas ^ 1- The distribution of stars, on 
the contrary, typically shows two peaks: one at ej* 1 and 
a second one at ej* ^ 0. This distribution is quite similar 
to tha t previously found by oth er authors (e.g. lAbadi et al.l 
l2003bl : IScannapieco et al.ll2009l l, and reveals the presence of 
two main components: a cold disc in rotational support and 
a hot spheroid supported by dispersion. 

Although the distribution of ej allows us to infer the 
relative importance of disc and spheroid, it is not enough 
for classifying individual stars into different components. In 
order to assign individual stellar particles to the different 
galactic components (thin and thick disc, bulge and halo), 
more parameters are needed in addition to the eccentricity 
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Figure 5. Left panel: Specific angular momentum at 2; = versus the mass for the dark matter haloes (open squares), baryonic discs 
(filled squares), and stellar bulges (starred symbols). Open circles are computed using the same methods than for observational galaxies: 
the specific angular momentum is obtained from the rotational curve, and masses from the mass-luminosity relation. Red symbols 
correspond to the low-resolution simulation. The limit radial distances for this plot were rvir for dark matter haloes, ropt = 3.2rs for 
discs and re for bulges. The solid and dashed boxes enclose the observational regions for spiral and elliptical galaxies, respectively, as given 
by Fall (1983). Right panel: Specific angular momentum versus rotation speed at 2; = 0. Small black dots correspond to late-type spirals 
compiled bv I Navarro! |1998). Filled circles correspond to the simulated galaxies, assuming that the rotation velocity is that measured at 
r = 2.2rs, and the angular momentum the sum of specific angular momentum of baryonic particles. Open circles indicate the location 
of the disc component, computed adopting the same angular momentum estimator as applied to observed data, j = 2rsVc. Starred open 
symbols represent angular momentum and rotation velocity of the stellar component. 



ej. In this work we will also use the total binding energy, 
computed using all particles inside the virial radius, and the 
ratio of the projected angular momentum of each baryonic 
particle on the plane of the disc, and the angular momentum 
expected if the same particle was in a circular orbit jp/Jc- 

As it can be seen in Fig. [71 where we plot all the stars 
in object LD-5003A, there exist heavily populated regions 
in the space defined by these three parameters. Those re- 
gions can be easily identified as the spheroid and the thin 
disc. In addition we also can recognize an intermediate, well 
populated region, which is also separated by means of these 
parameters and that we identify as the thick disc. 

In order to assign stars in different components we 
use an u nsupervised clus t ering algorithm, the ^kk- means ' 
method ISc holkopf etHI 1 19981 : iKaratzodou et al.1 l2004l : 
IPhillon et a l. 2004), which is the weighted kernel version 
of the 'k-means' algorithm. The 'kk-means' method uses the 
'kernel trick' in order to implicitly project all data into a 
non-linear feature space using a Gaussian kernel. Since the 
points are now located in a vector space of greater dimen- 
sionality, it is able, in principle, to capture clusters that are 
not linearly separable in input space, but are so in feature 
space. The algorithm attempts to minimize the sum of the 
distance in feature space from points to the assigned cluster 
centres, in the same way the standard k-means algorithm 
does. It proceeds iteratively from an initial selQ of cluster 

^ We run several instances of the algorithm to ensure it properly 
© 2011 RAS, MNRAS QQQ.fTti23l 
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Figure 7. Thick discs stars (green), thin disc stars (blue) and 
spheroid stars (red) decomposed by 'kk-means' in the 3D space 
defined by ej jp/Jc cind the total binding energy computed 
using all particles inside the virial radius, for simulated galaxy 
LD-5003A. 
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centres by assigning each point to the cluster with the closest 
centre. Then, the centres are re-computed as the centroids 
of each cluster. 

Given that the algorithm scales relatively poorly with 
the number of points and cluster centres, we do not ap- 
ply the algorithm to all stars in the main object. Instead 
we select a random sample of points to classify according 
to this method. The rest of stars are classified computing 
the Euclidean distance in the feature space, and selecting 
the nearest neighbours decided by majority vote. Finally we 
decide, studying the region spanned by each cluster, which 
group corresponds to the different components of the object. 
We have performed tests with different sample sizes to en- 
sure the sampling procedure does not affect the clustering 
results. 

The unsupervised nature of this segregation procedure 
implies that no prior, besides the expected number of com- 
ponents, have been assumed when separating the different 
components. In its construction we have also avoided the use 
of the properties (e.g., age or metallicity) expected for the 
stellar population of each component. Therefore, this decom- 
position method will allow for an unbiased study of many 
component properties and the connection with its history of 
mergers and accretion events in the cosmological context. 

In Fig. [8] we show the edge-on synthetic images of all 
stars in each galaxy, including satellites that will be excluded 
for the subsequent analysis, in top panels, and in subsequent 
rows we show the edge-on projection of stars in different 
components: thin disc, thick disc, and spheroid, respectively. 
The ej distributions of stars in each of these identified com- 
ponents are shown in Fig. [6l using the same colour code as 
in Fig. [71 In this figure, we can appreciate that the spheroid 
component (in red), shows almost no net rotation for most 
objects. Its distribution, however, tends to be skewed to- 
wards positive values of ej, this being specially clear for 
object HD-5101A. We interpret this skewness as a manifes- 
tation of a possible dual origin of the bulge, with a classical 
component supported by velocity dispersion, and hence a 
symmetric ej distribution, and a pseudo-bulge component 
with a net rotation (Obreja et al. submitted). In object HD- 
5103, A and B, it is clearly seen that stars and gas peak at 
different ej* values. This may indicate that the stellar disc 
is not aligned with the gas disc, probably due to the recent 
merger event in which it was involved. 

Table [3] summarises the main properties of each com- 
ponent. We see from this table that the fraction of stars 
assigned to the thick disc is ~ 16 per cent of the total 
stellar mass, in agreement with observational results (e.g. 
lYoachim DalcantonI l2006l : IColhns et al.l lioTol ). in which 
the fraction of stars locked in thick disc ranges from ^10 
to ^ 35 per cent. Columns (6) to (9) in this table, show 
the scale lengths of the thin (rsi) and the thick (rs2) discs 
obtained by fitting their r-band luminosity profile with a 
simple exponential, and the goodness of those fits (x^)- 
From these quantities we see that both disc components 
have similar radia l scale lengths, as in observed galaxies 
(|Yoachim DalcantonI l2006l ) . Columns (10) to (13) show 
the scale heights z^i and z^2 computed for thin and thick 



converges to a stable set of points that is independent of the initial 
choice of cluster centres. 



Table 4. Disc to total ratios measured from different methods. 
See the text for an explanation of the different methods used for 
the computation of these D/T values. 



Galaxy D/T? D/TPj^ D/TL 



LD-5003A 


0.70 


0.61 


0.34 


LD-5101A 


0.83 


0.75 


0.44 


HD-5004A 


0.68 


0.62 


0.39 


HD-5103A 


0.68 


0.56 


0.40 


HD-5004B 


0.76 


0.62 


0.36 


HD-5103B 


0.56 


0.57 


0.49 


HD-5004L 


0.72 


0.53 


0.42 



disc by just considering the fraction of RGB stars in a so- 
lar cylinder, as it will be explained in Section ISTTI We find 
that for low-density objects the mean values of thin and 
thick disc scale heights are lower than for objects in high- 
density regions, as we would expect. In any case, the ratios 
of the scale heights for the thick and thin discs in our simu- 
lated galaxies are found to be about 4 (excluding HD-5103), 
with a large spread, as usually f ound for external galaxies 
(e.g. lY oachim Dalcanton 2006), and i n our Galaxy (e.g . 
iGilmo re & Reid 1983; Siegel et al. 20^1; ljuric et al.l 120081 ). 
To end up, in the last column of Table [3] we display the asym- 
metric drift (Viag), computed as the difference between mean 
rotational velocity of thick and thin disc stars. These values 
show that the thick disc rotates at a slower pace than the 
thin disc as it occurs in our Galaxy (Soubiran et al. 2003r) 
and in some external galaxies (|Yoachim Dalcantonll2005l l . 

Finally, in Tabled we show different D/T ratios mea- 
sured with different techniques. D/T^ represents values ob- 
tained from ID fits of the total i-band luminosity profile. 
The next column shows the D/T^j^ ratios obtained when the 
luminosity profile of particles in the 'kk-means' spheroidal 
component is fitted with a Sersic function and, separately, 
the luminosity profile of particles in the 'kk-means' disc com- 
ponents is fitted with a single exponential. Finally, D/T^j^ 
gives the direct ratios, without any fit, of the 'kk-means' 
disc to total luminosity. We see from this table that the two 
fitted ratios (D/T^ and D/T^j^) have roughly similar val- 
ues for each object, with a systematic tendency for D/T^k 
to be slightly lower than D/T^. On the contrary, the non- 
fitted direct ratios, D/T^k, of the 'kk-means' components are 
significantly lower than D /T^. These results are i n agree- 
ment with those found in Scannapie co et al ] (|2010l ). where 
our fitted and non- fitted D/T ratios must be compared to 
their photometric and kinematic D/T ratios, respectively. As 
noted by those authors, the difference obtained from these 
two methods can be partly explained from the fact that, in 
the 'kk-means' (or kinematic) decomposition, disc stars are 
not present or hardly present within the inner regions, while 
the fitted photometric decomposition assumes an exponen- 
tial profile starting at the very centre of the galaxy. These 
results show that, although the photometric D/T ratios can 
be employed to compare the properties of simulated and real 
galaxies, kinematic properties of galaxies are not prop erly 
inferred from such D/T ratios (jScannapieco et aLll201Ql V 
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Figure 8. Edge-on maps of stellar magnitudes (r-band) for the simulated galaxies at z = 0. All images are 40 kpc side. Top panels show 
the total stellar content of each object, including satellites that will be excluded for the analysis, and in next rows we draw, from second 
row to bottom, thin stellar disc, thick stellar disc and spheroidal component as classified by our method. 



Table 3. Main properties of components: Stellar mass enclosed in the spheroid (Mgph), thin disc (Mthin), and thick disc (Mthick)- 
Fraction of mass in the thick disc compared to the total stellar mass in the object (fthick)- Scale lengths of the thin (r^i) and thick (rs2) 
disc components, and the goodness of the fits (x^). Scale heights of the thin (z^i) and thick (2^2) discs, and the goodness of the fits 
(x^). Finally, the asymmetric drift (Viag). 
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(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


LD-5003A 


1.00 


3.42 


1.95 


0.13 


2.95 


0.0022 


2.87 


0.0001 


0.32 


0.0291 


1.49 


0.0106 


-12.31 


LD-5101A 


0.52 


3.17 


1.69 


0.17 


3.83 


0.0032 


3.45 


0.0002 


0.45 


0.0017 


1.35 


0.0102 


-4.87 


HD-5004A 


1.86 


5.98 


5.43 


0.18 


4.16 


0.0011 


3.81 


0.0001 


0.42 


0.0108 


2.33 


0.0015 


-16.53 


HD-5103A 


1.73 


5.65 


5.42 


0.19 


4.19 


0.0016 


4.26 


0.0006 


0.67 


0.0061 


8.32 


0.0004 


-27.75 


HD-5004B 


1.74 


5.79 


3.89 


0.14 


3.58 


0.0027 


3.39 


0.0014 


0.34 


0.0059 


2.17 


0.0050 


-5.41 


HD-5103B 


1.65 


7.80 


4.65 


0.16 


3.81 


0.0002 


4.76 


0.0004 


1.01 


0.0005 


13.4 


0.0005 


-28.51 



4.2 Observational-like samples 

A full decomposition method like that described in the pre- 
ceding Section 14. 1[ where each stellar particle is assigned 
to a galactic component, is very useful for the modelling of 
the galactic structure and formation processes, as well as 
for comparison with other simulation-based studies. How- 
ever, it can not be applied to observational data. In order 
to facilitate the comparison of our results with those avail- 
able from observations, we have also selected spheroidal and 
disc subsamples based on two different criteria similar to 
those frequently used in observational studies. Each of these 



observational- like samples are not intended to contain all 
the stellar particles of a given component, but just a subset 
of particles easily identified in observations and that mostly 
belong to such a galactic component. 

The first criterion is mainly founded on the positions of 
star particles in the galactic object. The bulge sample con- 
sists of particles in the innermost region (r ^ 5re) above 
and below the galactic plane {2zhi ^ \z\ ^ ^Zhi) to min- 
imize the contribution of disc stars. The two disc samples 
are selected in a cylindrical shell 2rs ^ r ^ Sr^, which min- 
imizes the contribution from bulges, and imposing a min- 
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imum rotation speed of 50 km s ^ in ord er to minimize 
the contribution from halo stars (jSales et al. 2009). The dis- 
tance to the galactic plane determines the difference of both 
samples: stars located close to this plane {\z\ ^ 2zhi) con- 
stitute the thin disc sample, while those at farther heights 
{zh2 ^\z\ ^ Szh2) constitute the thick disc sample. 

The second method uses kinematic information to select 
disc stellar samples, as employed in Galac tic stellar studies 
(jSoubiran et al.(l2003l : iBensbv et alll2003l ). In those works, 
thin and thick disc stars are selected by assuming that the 
Galactic space velocities of such stellar populations have 
Gaussian distributions 



f{U,V,W)<xexpi--^ 



{V-{V)f 



2al 



(3) 



where U, V and W are the radial, azimuthal, and axial com- 
ponents of the stellar velocity, respectively, and au , cfv and 
aw the velocity disp ersions found for thin and thick disc 
Galactic stars (see e.g. lSoubiran et al]l20 0j). Knowing these 
values, stars in the Galactic disc are classified according to 
their major probability of belonging either to the thin or to 
the thick disc. In our simulations, however, we can not use 
the values of the Milky Way velocity dispersions, being nec- 
essary to compute them. To do that, we first selected a sam- 
ple of stars in 'solar neighbourhood volumes', equivalent to 
cylindrical shells normal to the disc plane between two and 
three scale lengths of the thin disc (2rs ^ r ^ 3rs), which 
will be mainly populated by disc stars, and then we segre- 
gated them in two components according to their kinematics 
and metallicity, as in Soubiran et al. (2003). These two new 
samples are employed to compute the characteristic velocity 
dispersions that we will apply in the classification of stars 
in the disc, not only in th e ^ sola r neighbourhood volume'. 
According to lBensbv et al ](l2003'), in order to minimize the 
contribution of the thin disc stars in the thick sample (and 
vice- versa), we classify as thick disc stars those which are 
at least ten times more likely to be a thick disc star than a 
thin disc star, i.e., those with frn/fD ^ 10. Similarly, thin 
disc stars are selected as those with fTD/fo ^ 0.01. 



5 VERTICAL DENSITY PROFILES 

5.1 Scale heights of the disc components 

The variation of stellar density with the distance z, above 
or below the Galactic plane, can be determined by count- 
ing faint stars in fields centred on the Galactic poles. Stellar 
counts in the solar neighbourhood result in vertical profiles 
which can be fitted wit h a double expon e ntial function. This 
fact was interpreted bv lGilmore k ReidI ('1983) as due to the 
presence of two different disc components, thin and thick, 
embedded in the Galaxy. Values measured for this double 
fit rises to ^ 350 pc for thin disc scale he ight, and ranges 
from 750 to 1450 pc for the thick one (e .g iGilmore ReidI 
ll983l : ISiegel et al.ll2002l : ljuric et al.l l2008'). In external galax- 
ies, however, resolving stars is only possible for a small num- 
ber of galaxies, but the presence of thick discs have been pho- 
tometr ically detected as an excess flux at large galactic lat- 
itude s (jPalcanton &: Bernstein|[2002l : lYoachim &: DalcantonI 
12006). These works show that thick discs are ubiquitous in 
disc galaxies with scale heights spanning a wide range of val- 
ues, but larger than the scale height of the embedded thin 



discs by a factor 2. llbata et al ] ("2009') measured the scale 
height for the Milky Way analogue NGC 891 from stellar 
counts and obtained a value, Zh = 1-44 ± 0.33 kpc, quite 
similar to that observed in the Galaxy. 

In order to compute the vertical density profiles of simu- 
lated galaxies, we have followed a procedure similar to that 
commonly used in observational studies. The contribution 
of spheroidal components was minimized by just consider- 
ing stellar particles with projected galactocentric distances 
in the range 2 < r/rs < 3. We then counted the number 
of RGB stars implied, for the assumed IMF, by the age 
and metallicity of the stellar particles lying within differ- 
ent vertical bins. The vertical profile of the surface density 
was then obtained for each simulated galaxy by averaging 
several edge-on projections. Such profiles are shown (black 
solid lines) in the left panels of Fig. [9l where we also rep- 
resent the contribution of the thin (blue squares) and thick 
(green triangles) discs and the spheroid (red pentagons) de- 
composed by 'kk- means'. The corresponding exponential fits 
of both disc components are also shown (red lines) in Fig. [9l 
and their values are summarised in columns (10) and (12) 
of Table (3] Fits of the thin disc take into account only the 
innermost part of the profile, similarly the thick disc fits 
ignore the innermost part, dominated by the thin disc. 

The scale heights obtained in our low-density objects 
have mean values of ^ 0.4 kpc0 for the thin disc and ^1.4 
kpc for the thick component. In the high-density objects 
these scales are somewhat higher, specially those of the thick 
discs, which have scale heights significantly larger than 2 
kpc. This is likely due to the more active history of mergers 
and accretions for galaxies in crowded regions. Indeed, the 
object HD-5103A (and its counterpart HD-5103B with less 
SF efficiency) has suffered a recent merger and for this rea- 
son it is perturbed at z = 0. As can be seen in the figure its 
vertical profile is almost exponential and it is dominated by 
thick disc stars, even at farther heights where halo stars are 
supposed to be the dominant population (red pentagons in 
Fig. [9|). In any case, the ratios of the scale heights for the 
thick and thin discs {zh2/zhi) in our simulated galaxies are 
found to be about 4, with a large s pread, as usually found 
for e xternal early spirals (see, e.g., lYoachim &; DalcantonI 
l2006l . focused on Sd galaxies, but containing a compilation 
for SO-Sb galaxies). 

5.2 Radial variation of the disc scale heights 

Another important feature of observed thick discs is that 
their vertical profiles do not significantly depend on the 
galactocentric distance, i. e., the scale height seems to be 
roughly constant with radius (Ro bin et al ]"2003";' Iba ta et al.l 
2009 ). Such a result is difficult to be reconciled with a 
merger-induced scenario for thick di sc formation, which in - 
evitably leads to fiared thick discs fBournaud et al.l l2009l V 
and favours models where the disc thickening results from 
internal processes such as gravitational instabilities. 

In order to test whether the scale height changes 
with radius , we h ave followed the procedure employed by 
llbata et al.1 (|2009l ) for an observed external galaxy. We have 

^ Note that thin disc scale height is close to our gravitational 
resolution, and thus it should be regarded with caution. 
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Figure 9. Vertical density profiles of all stars within two vertical cuts of same width and different distance to the galactic centre. Black 
solid lines, blue squares, green triangles and red pentagons represent the total, thin disc, thick disc and spheroid stellar counts in each 
section, respectively. Red solid lines at the left panels stand for the fits of thin and thick discs in the outermost region. Red dashed lines 
at the right panels show an exponential fit to the profiles in the internal cut constrained to have the same scale heights for both, the thin 
and thick disc, as those respectively measured in the outermost region. 



drawn (right panels of Fig. [9]) the vertical density profile of 
stars within a region (1 < r /vs < 2) inner than in the pre- 
ceding analysis. For this internal region we have computed 
the exponential fits of both discs components constrained to 
have the same scale heights, zui and Zh2^ as in the external 
one. We find that, for the simulated objects (LD-5003A and 
LD-5101A) in low-density regions of the Universe, roughly 
constant Zh values would result in acceptable fittings for 
both the inner thin and thick. By contrast, for objects in 
high density-regions (HD-5004A and HD-5103A, and their 
counterparts HD-5004B and HD-5103B with less SF effi- 
ciency) , roughly constant Zh values give acceptable fits only 
for the thin disc of some simulated galaxies, but clearly fails 
for the thick disc. Therefore, for these latter objects we can 
not reject the merger-induced scenario as a key ingredient 
for thick disc formation. 



6 STELLAR AGES 

6.1 Age distributions of the different components 

Fig. [10] shows the distribution of stellar ages for the differ- 
ent components of our simulated galaxies. Each distribution 
was normalised to the total stellar mass of the correspond- 
ing individual component. The different rows in this figure 
depict the age distribution of samples selected from the dif- 
ferent criteria described in Section O the first row refers to 
'kk-means' decompositions, the second row to observational- 
like samples built from position criteria, and the third row 
to observational-like samples built from kinematic criteria. 
The corresponding median values of stellar ages are sum- 
marised in Table O while the star formation histories for 
each object and its 'kk-means' components are shown in 
Fig. [TTl Note that, for runs LD-5101A and HD-5103B, the 
observational- like samples for the bulge component were too 
small for statistics and, hence, they are not depicted in Fig. 

[ini 

A clear trend is found in all simulated galaxies for any 
kind of sample. The thin disc is always the youngest com- 
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Age(Gyr) 

Figure 10. Stellar age distribution normalised to the total mass of each individual component. In each row we depict the age distributions 
obtained with different classification methods, as indicated in the y-axis. In the first row we show age distributions for star particles in 
each component identified by our classification method. In the second one samples are selected by positions in the simulated galaxies, 
and in the last one, kinematic criteria have been employed to identify stars in thin and thick disc. In all panels, blue line refers to thin 
disc star particles, green is for thick disc and red line stands for the spheroid. 



ponent with median ages ranging from 3.8 to 6.7 Gyr. Just 
a small fraction of the thin disc stars have ages older than 
10 Gyr, while just a small fraction of stars in the spheroid 
are younger than 5 Gyr. The thick disc is instead populated 
by stars with intermediate ages. These results are consistent 
with the observational estimates for the different compo- 
nent s in the Milky Way and nearby galaxies. Indeed, several 
f e. g. iFuhrmannll 1 998l : Bensbv et al ] |2003l . l2005l : lReddv et al.l 
agree that stars in the thick disc of the Galaxy have 
mean ages greater than 8 Gyr meanwhile the estima tes for 
the thin disc stars lay between 5 to 7 Gyr. Bensbv ~t al.l 
(|2QQ5l ) derived a mean age of 4.3 ± 2.6 Gyr for stars in 
the thin disc of the Galaxy, and 9.7 ± 3.1 Gyr for a sam- 
ple of thick disc stars. Similarly, for nearby galaxies, pho- 
tometric studies find that the off-plane regions are domi- 
nated by old stars (> 6 Gyr) (^Yoachim Dalcanton.2006t : 
iDalcanton k, BernsteinI l2002l l . In their study of the inte - 
grated spectra of six galaxies, lYoachim DalcantonI (|2008l l 



obtained a median age of 4.6 Gyr for the thin disc, and 
ages ranging from 3.8 to 10.9 Gyr ( with a median va lue of 
7.2 Gyr) for the thick disc. Scanna pieco et al.l (|2011bl ) have 
also found similar results for the study of stellar populations 
in their simulated galaxies. The agreement between results 
obtained with different simulation techniques and decompo- 
sition methods, could be interpreted as a general feature in 
the age of stellar populations of simulated galaxies. 

Focusing on star formation histories, we can see in Fig. 
[11] that the thin discs of all simulated galaxies are currently 
forming stars with a star formation rate lower than about 
one solar mass per year. For almost all simulations the bulk 
of stellar mass formed 10 Gyr before the present epoch, pop- 
ulating preferably the spheroid. At low redshift, however, 
there is just some residual star formation in this component. 
On the contrary, the thick disc is not forming new stars at 
low redshift s. The stellar particles settled in this component 
were formed during an extended epoch in the past, with a 
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Table 5. Stellar median ages in Gyr for different components of 
each simulated galaxy. The three rows for each object correspond 
to a 'kk-means' decomposition, samples from position criteria, 
and samples from kinematic criteria, respectively. 



Galaxy 




Ap^ej-u- 1 
-'^fe'^tnick 




LD-5003A 


9.6 ±2.2 


8.7 ±2.3 


4.7 ±3.3 




10.7 ± 1.0 


9.0 ±2.0 


4.3 ±3.1 






8.5 ±2.3 


4.1 ±2.9 


LD-5101A 


10.8 ± 1.9 


8.1 ±3.0 


3.8 ±3.3 






8.1 ± 2.4 


3.3 ± 3.2 






7.2 ±2.9 


3.1 ±3.0 


HD-5004A 


9.7 ±2.2 


8.3 ±2.5 


4.7 ±3.3 




10.3 ± 1.6 


8.3 ±2.4 


5.0 ±3.6 






9.5 ±2.5 


5.3 ±3.4 


HD-5103A 


9.0 ±2.9 


9.0 ±2.8 


5.7 ±3.8 




9.4 ±2.7 


9.8 ±2.3 


3.1 ±3.8 






8.6 ±2.9 


3.8 ±3.6 


HD-5004B 


9.5 ±2.2 


8.6 ±1.9 


4.9 ±2.7 




10.2 ± 1.5 


9.3 ±1.8 


4.3 ±2.7 






9.1 ±2.2 


4.3 ±2.6 


HD-5103B 


8.8 ±2.0 


8.7 ±2.9 


6.7 ±3.7 




9.9 ±1.8 


9.4 ±2.7 


6.2 ±3.8 






9.2 ±2.8 


6.4 ±3.7 



low rate similar to that observed in the thin disc. In the 
last 10 Gyr merger and accretion events triggered the star 
formation leaving a trace as peaks in the star formation his- 
tory. This is clearly visible in simulation HD-5103A and its 
counterpart HD-5103B, which have been involved in a re- 
cent major merger. The new stars formed during this event 
mainly increase the stellar content in the spheroid and the 
thin disc. Differences in the SF parameters produce differ- 
ences in the response of simulated galaxies to merger events. 
In those objects in which SF has been difficult, the SF oc- 
curs more quietly and during a more extended epoch than in 
their counterparts with a more efficient SF, in which those 
events produce sharper peaks. 



6.2 Age and colour gradients 

Colour profiles are related with galactic properties such 
as metallicity, stellar age or dust content, and for this 
reason they have long been used for exploring galactic 
properties in large galaxy samples and catalogues (e.g. 

120101 1. Fig. [T2l shows the g-r colour 



iGonzalez-Perez et al 



and age profiles for each simulated galaxy. Black solid lines 
depict the total profiles, while coloured dashed lines show 
the individual profiles of each galactic component identi- 
fied by the 'kk-means' method. From this figure we see 
that the colour profiles are dominated by the thin disc 
contribution, which is responsible for their U shape: every 
galaxy and thin disc in our sample is redder in the inner- 
most region and becomes bluer as we move away in the ra- 
dial direction until a minimal value is reached. Then, the 
stellar content becomes progressively redder. These results 
are in good agre ement with obser v ations of Type-II spi- 
ral galaxies (e.g. Bakos et al. l2008l : lYoachim DalcantonI 
120081 : iGonzalez-Perez et al.ll2010l ) , where such U-shaped pro- 
files are also found, with median colour gradients close to 
zero. These features have been recen tly studied using cosmo- 
logical hydrodynamic simulations ([Mart mez- Serrano et al.l 




5 10 5 10 5 10 
Lookback Time (Gyr) 

Figure 11. Star formation history for each component. Black 
solid line stands for the total star formation rate in the object. 
Thin disc is represented with a blue dotted line, and green short 
dashed line for the thick disc. Star formation rate in the spheroidal 
component is plotted with a red long dashed line. 



I2OO9I : ISanchez-BlazQuez et al.l l2009l ) and tend to favour an 
inside-out formation scenario for the galactic discs. Accord- 
ing to such a scenario, star formation has moved progres- 
sively towards outer regions of galactic discs. As a conse- 
quence, older and more metal-rich stellar populations are 
located in the innermost regions of discs leading negative 
colour gradients. Due to the secular evolution, stars mi- 
grate from the inner to the outer parts of the galaxy, 
reaching the outskirts where the specific star formation 
is very small. These outskirts progressiv ely become pop- 
ulated by old and red stars (Mar tinez- Serrano et al. I I2OO9I : 
ISanchez-Blazquez et al' 

In contrast to that found for the thin disc, the colour 
and age profiles depicted in Fig. [12] for the thick disc and 
spheroidal components are quite uniform across the galaxy, 
with little sign of radial gradient. These profiles are a clear 
indication of a different formation scenario for such com- 
ponents: their stellar content is older and was formed in a 
shorter time than in the thin disc (see also Fig. [10] of age 
distributions). As in observed galaxies (MacArthur et all 
age gradients in the spheroid are almost zero, with 
just small deviations from a flat profile both towards nega- 
tive or positive values. These deviations are probably more 
visible for objects in low density regions, where the spheroid 
shows a slightly negative age gradient. 



7 CHEMICAL PROPERTIES 

7.1 Metallicity Distribution Functions 

Element abundances of individual stars and gas provide a 
valuable record of the star formation and enrichment history 
during the early stages of galaxy formation and evolution. 
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Figure 12. Top panels: Colour profiles for each component in the galaxy as seen face on. As in previous figures solid line stands for 
the colours of all stars in the galaxy. Thin disc stars are plotted using a blue dotted line, and the colour gradient of thick disc stars is 
shown with a green short dashed line. Red long dash ed line stands for c olour gradients in spheroid stars. For reference, the observed 
colour profiles of 39 Type-II galaxies, as reported by iBakos et al.l (|2Q08l ) are shown as a shaded grey area. Note that the population 
synthesis used in this paper does not include dust effects. Bottom panels: Mean stellar ages weighted by r-band luminosity for each 
simulated gal axy. Line code is as in top panel and previous figure. Similar profiles weighted by mass can be found in our previous paper 
([Martinez-Serrano et al.ll20Q9l ). 



and it has long been used to propose or tes t the different 
me chanisms involved in such processes (e.g. iGilmore et al.l 
[1995: Reddv et al. 20061: ISchuster et al.|[2006l ). The fact that 
kinematic and chemical properties of thick disc stars differ 
form those of thin disc stars has been interpreted as an ev- 
idence of their distinct nature, making possible the differ- 
entiation of stars from both components using their chemi- 
cal properties (e.g. Schuster et al. 2006; Schonrich & Binne^ 
[2009bh . Otherwise iMelendez et al.l (|2008l ) found no signifi- 
cant abundance differences between the bulge and the lo- 
cal thick disc stars in their study, which is interpreted as 
a consequence of the rapid star formation history of both 
components, rather than any deeper relation between them. 

In the Milky Way, thick disc stars have been found to 
be old and more metal-poor and_also_a- enhanced co mpared 
to stars in thin disc (e.g. Fuhrmann lIS^S ; Bens bv et al.l 
l2003l . l2005l : iReddv et al.ll2006l : iRuchti et al.ll2010l ). To con- 



strain the chemical evolution of the Galaxy, samples of long- 
lived stars, as F an d G dwarfs, have been employed (e.g . 
Gilmore et all Il995l: iFuhrmannI Il998l : ISensbv et all 12001 
20051 : iReddv et al.ll2006l ). The metal distribution of thick 
disc stars is roughly characterized as a Gaussian with mean - 
iron abundance ^ -0 .7 to -0.5 dex (e.g. iGilmore et al.|[l995l : 
ISchuster etHI 120061 ). with most stars between -0.5 to - 
1.0 dex, but with a long tail of metal-poor stars which 
extends to about -2.2 dex, values usually associated with 
halo stars. This metal-poor tail is interpreted by some au- 
thors as an evi dence of the complex structure within th e 
thick disc (e.g. 'Schuster et al.l l2006l : ICarollo et aD I2OIOI I. 
Yoachim & Dalcanton (2008) used Lick metallicity-sensitive 
indices to derive metallicity properties of thin and thick disc 
components in six external galaxies. They found metal-poor 
and no significant a-enhanced thick discs, probably due to 
the insensitivity of Lick indices to the differences found at 
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low metallicities. The median metallicity value measured in 
the thick and thin disc of these galaxies is [Z/H]=-0.6, as 
expected for low-mass galaxies. 

In Fig. [13] we depict the MDF of each component in our 
simulated galaxiefl. As in Fig. [10] each row corresponds to 
samples selected from the different criteria described in Sec- 
tion 3] and each distribution is normalised to the total stellar 
mass of the corresponding individual component. In general 
we find that the thin disc is the most metal-rich component, 
with median values close to the solar value, between -0.19 
to -0.06. Moreover, the distribution is quite similar for the 
different thin disc samples. The spheroid, on the contrary, is 
the most metal-poor distribution, with median values cover- 
ing a wide range between -0.61 to -0.27. In this component 
a clear difference is observed among the different samples. 
In those based on the 'kk-means' method we find a bimodal 
MDF, while such a bimodality is not observed in the sam- 
ples selected from position criteria, which exclude stars near 
the galactic plane. Since the distribution of parameters ej 
and jp for these high-metallicity stars in the spheroid clearly 
corresponds to this component, such a high-metallicity peak 
seems to be a real feature of the spheroid MDF, and not an 
artefact due to disc stars that eventually could be misclassi- 
fied as belonging to the spheroidal component. The bimodal 
MDF for the spheroid seems more likely due to the differ- 
ent formation mechanisms of the bulge of spiral galaxies: 
the low- metallicity peak would correspond to the classical 
bulge, while the high-metallicity peak would correspond to 
the bulge fraction formed from instabilities of the inner disc 
(Obreja et al. submitted). 

With the only exception of simulation HD-5103A and 
its counterpart HD-5103B, in which the stellar content is 
well mixed, the thick disc MDF shows an intermediate be- 
haviour between the thin disc and the spheroid. This compo- 
nent has in fact a metal content poorer than in the thin disc, 
while richer than in the spheroid. Their MDFs peak at val- 
ues from -0.45 to -0.24, but their shape is not as sharp as for 
the thin disc and extends to very metal-poor regions. In the 
case of observational- like samples selected from kinematic 
criteria (third row in Fig. [T3|) the thick disc MDF shows a 
clear bimodality, which could be relate d with the metal-poor 
thick disc defined by some authors ([Schuster et al.1 l2006l : 
ICarollo et al]l2010l l. 



7.2 Metallicity variation along the disc 

The distribution of element abundances within galaxies is 
another important feature in the comprehension of galaxy 
formation and evolution. First measures, based in emission 
lines in large samples of HI I regions in spiral galaxies (e.g . 
Shaver et al] 'l983: 'Zaritsk v et al1 ll994': 'Garnett et allll997l : 
van Zee et al . 1998), found radial abundance gradients in al- 
most all large spi ral galaxies, with in ner HI I regions richer 
than outer ones. I Shaver et al.l (|l983l l measured an oxygen 
abundance gradient of ~ —0.07 dex kpc~^ in the Milky 
Way, a value quite similar to that obtained bv I Luck et al.l 
(|2006l ) fitting iron abundances of young stars in the Galaxy, 



^ The solar abundances of lCrevesse Sauvall (|l998l V with Zq = 
0.0169 have been assumed for the subsequent analysis. 



or lPancino et al 1 (|2010h using open clusters. This result sug- 
gests that the chemical evolution of the Galactic disc is de- 
termined mainly by local chemical evolution, because merg- 
ers would have produced more complex abundance distri- 
butions. In external regions, however, steeper gradients are 
found, wit h metal-poor stars at high radial distances (e.g. 
lYong et"al. 2005). The existence of such a radial profile can 
indicate that the unobserved inner stellar thin disc may be 
significantly more metallic than the solar neighbourhood. 

Observations of the Milky Way reveal that the differ- 
ent g alactic components ha ve characteristic abundance pat- 
terns. iNavarro et al.l (|201ll ) show distinctive abundances for 
the thin and the thick disc in the solar neighbourhood, with 
the thick disc being significantly less enriched. The spheroid 
is known to have a gradient with the outer (halo) parts being 
very metal poor (jlvezic et al.ll2008l ) , and the inner bulge hay - 
ing metallicities around the solar value (|Zoccali et al.ll20'08l V 

In Fig. [14] we depict the metallicity vs. radius distri- 
butions for each component identified by the 'kk-means' 
method in our simulated galaxies. In addition, the first 
row represents the gas component, with computed gradients 
ranging from —0.060 dex kpc~^ in the object HD-5004A to 
—0.039 dex kpc~^ in LD-5101A which is the less enriched 
galaxy. 

At first glance, one can appreciate that the more 'disky' 
a component is, the lower its dispersion in metallicity for any 
given radius, with the spheroids showing the highest range 
of variation in metallicity, and the thin discs the lowest. The 
thick discs appear here as an intermediate component, with 
almost the same radial extent as the thin disc but a higher 
dispersion in their metallicity values. The least massive ob- 
jects (LD-5003A and LD-5101A) show the lowest metallicity 
gradient for the thin disc component, whereas a more mas- 
sive galaxy (HD-5004A and its counterpart HD-5004B) show 
a larger gradient, with the run having the least effective set 
of SF parameters (HD-5004B) having a larger gradient. 

Given that our sample is very small, no conclusions can 
be drawn with respect to the relation between morphological 
type (or D/T ratio) and metallicity gradie nt, although there 
is a s imilar effect known from observations ([Edmunds Rovl 
I1993I I. The difference between HD-5004A and HD-5004B, 
given that it is the same object with a different set of SF 
parameters, is however more telling. A reduction in the over- 
all SF efficiency leads to a steeper gradient. Given that the 
metallicity values in the inner part of both objects is sim- 
ilar, the difference must be due to a reduction in the SF 
efficiency in the outer parts of the disc, where the density 
threshold criterion is not always met by the HD-5004B simu- 
lation. This effect is present in obser vations of nearby spirals 
(|Lerov et al.ll2008l :lBigiel et al.l l201Qh . where a decreasing SF 
efficiency with density (or galactocentric radius) is found, 
leading to extremely inefficient SF in the outskirts of discs. 

The spheroids show a rather complex structure in Fig. 
1141 For most objects several bumps in the MDF of the bulge 
are apparent (see also Fig. [13]). As was explained in the pre- 
ceding subsection, this is likely the outcome of two different 
processes driving the formation of the bulge, with the high- 
est metallicity peak being associated with a pseudo-bulge 
and the lower metallicity distribution being associated with 
the classical bulge. Further results on the bulge abundances 
and their formation process will be presented in Obreja et 
al. (submitted). 
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Figure 13. Metallicity distribution function of stars in each galaxy component. As in Fig. llOl in each row is depicted the MDF of samples 
selected from the different described criteria. Each component is represented by different colours, as in previous figures. 



Object HD-5103A, and its counterpart HD-5103B, de- 
serve an special mention. Since it has undergone a recent 
merger, its stellar content is mixed, and the distributions are 
generally broader for any given radius, with many particles 
shifting orbit and therefore lower particle densities attained 
in the metallicity vs. radius plane. 



7.3 a-elements 

A number of works in the literature, using different stel- 
lar samples, have studied the abundance of a-elements in 
the bulge, thin dis c and thic k disc of the Milky Way (e.g. 
Bensbv et alll2003l. [20 05: Red dv et al.ll2006l: Ijablonka et al.l 
20071: iMelendez et al.l l2008>: .Ruchti et al.l l2010ll. Accordine: 
to these works, the Galactic bulge and thick disc are on 
average more enhanced in a-elements than the thin disc, 
suggesting a short formation time-scale for these two com- 
ponents. The a-elements are in fact mainly produced during 
Type II SNe on a short time-scale, whereas iron is also pro- 
duced by Type la SNe on a much longer time-scale. An en- 
hanced abundance in a-elements can be then interpreted as 



a signature of a short period of stellar formation, where SNe 
la have no enough time to dilute the a-elements produced 
by Type II super novae. 



Moreover, in the thin and thick disc the abundance of 
a-elements show tight and distinct trends with metallicity. 
In the thin disc, this trend has a gentle slope wherea s, for 
the thick disc, some authors (jBensbv et al.|[2003lj2005l l have 
pointed out the presence of a down-turn (or 'knee'). Accord- 
ing to these authors, for metallicity values [Fe/H]< —0.4, 
stars in the thick disc have [a/Fe] values higher than in the 
thin disc. At [Fe/H]^ —0.4, the [a/Fe] trend has a down- 
turn, so that it declines toward solar values merging with the 
thin disc trend for [Fe /H ] > —0.4. In contrast, other authors 
(e.g. Red dv et all 12006') doubt the existence of this 'knee', 
but agree in the different behaviour of thin and thick disc 
stars. 



In Fig. [15] we show the distribution of the abundance of 
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Figure 14. Metallicity radial distributions for gas, thin disc, thick disc and spheroid (from top to bottom). 



a-elements represented by oxygerQ. In general, all samples 
of thick disc stars distribute around a median value 0.3) 
slightly higher than in the thin disc. For the location-based 
samples of stars in the spheroid, we also find a-enhanced 
abundances, with median [0/Fe] values (~ 0.45) higher than 
in the thick disc. The spheroid samples identified by 'kk- 
means' have instead a bimodal distribution: there exists a 
population of stars poor in a-elements, as well as another 
population enhanced in these elements. Such a bimodality 
could be due to the combined action of two bulge formation 
mechanisms (classical bulge and pseudo-bulge), and also to 
diff erent origin of stars in the halo as recently pointed out by 
iNisse n & Schuster (2010). They suggest that low-a stars in 
the Galactic halo probably were accreted from dwarf galax- 
ies that had lower star formation rates, while the high-a 
population is formed by the first stars likely originated in 
a dissipative collapse of the proto- Galactic gas cloud, or by 
ancient disc or bulge stars heated to the halo by merging 
processes. 

In Fig. [16] we depict [0/Fe] versus [Fe/H] for disc 
stars in each simulated galaxy. To com pare our resul t s with 
the o bservational trends reported by Bensbv e t al.l (|2003l , 
I2OO5 IV the number of points in Fig. [16] are proportional to 
the fraction of F and G stars in each star particle, com- 
puted from its age and metallicity assuming a Salpeter 
IMF. From this figure we can see that at low metallici- 
ties the thick disc stars are generally more a-enhanced than 
those with the same metallicity in the thin disc. This is 
particularly clear in objects LD-5003A and LD-5101A, lo- 
cated in sparsely populated regions of the universe. In these 
objects, a 'knee' is found at roughly the same metallicity 
value ([Fe/H|^ ~^-^] P Milky Way observations of 
iBensbv et all (|200l I2OO5I ). We also find for these objects 



^ Similar trends are found when the abundance of a-elements is 
represented by magnesium 



that there are some old thin disc stars over-abundant in 
oxygen. Although less clearly, such an a-enhancement of the 
thick disc low-metallicity stars seems to be also present in 
the objects HD-5004A, HD-5103A and their counterparts, 
populating denser environments. We also note that, except 
for HD-5103B, the plot of [0/Fe] versus [Fe/H] has sec- 
ondary branches containing stars of both disc components 
with high a-abundance and metallicity. These secondary 
curves are mainly populated by young thin disc stars formed 
in an enriched medium and by old thick disc stars which 
formed during a rapid episode. These old stars enriched the 
medium with a-elements and then they have been evolving 
at the same time that the iron abundance has been increas- 
ing. 



8 SUMMARY AND CONCLUSIONS 

We have presented a detailed analysis of the photometric 
and chemical properties of seven simulations of disc galaxies 
in the ACDM scenario. Each simulation is a cosmological 
zoom-in that includes high-resolution gas and dark matter 
particles for the converging region that generates the main 
object. Star formation is modelled through a Kennincutt- 
Schmidt-like law with a star formation efficiency that implic- 
itly accounts for energy feedback. Although the addition of 
some explicit form of kinetic feedback could presumably lead 
to simulated gala xy populations in even better agreement 
with observations fceres et al.ll2009l \ in this paper we have 
preferred to test the minimal conditions for the formation 
of realistic discs, thereby facilitating the physical interpreta- 
tion of our results. The methods for chemical feedback and 
abundance-dependent cooling are also included. The simu- 
lations presented in this paper correspond to four different 
objects selected from both sparsely populated and crowded 
regions. In order to test the robustness of our results, three 
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Figure 15. Over-abundance distribution of stars in each galaxy component measured as the oxygen ratio over iron. 



of the seven runs consisted of variants with different SF ef- 
ficiency or resolution. 

From the bulge-disc decomposition of their luminosity 
profiles, we find that our simulated galaxies are disc domi- 
nated, with bulge sizes comparable to those observed in late- 
type spirals. In addition they have also realistic kinematic 
properties and lie in the region defined by the observational 
Tully-Fisher relation. 

In order to analyse the fine structure of these objects, 
we have performed a dynamical decomposition of the galac- 
tic components based on the total binding energy and the 
axial and radial components of the angular momentum, nor- 
malised with the angular momentum expected for a particle 
in a circular orbit. We used an unsupervised clustering al- 
gorithm which classifies stars by minimising the sum of the 
distances of each particle to its cluster centre in the space 
defined by these three parameters. In this way, we can iden- 
tify stars from the thin and thick disc, and also from the 
spheroid, which includes bulge and halo. We will now dis- 
cuss our main results for each component: spheroid, thin 
disc and thick disc. 



We find that the spheroidal component of the simu- 
lated galaxies is composed mainly by old stars, with just a 
small contribution of stars younger than 5 Gyr. The bulk of 
stars formed in a rapid collapse at very early times. After- 
wards, the star formation turned off, and only residual star 
formation is observed at low redshifts. The age and colour 
profiles for this component show that the stellar content 
is distributed quite uniformly, as observed in real galaxies, 
with mean ages almost constant along the galaxy. The only 
exception is the simulated object HD-5103A and its less effi- 
cient SF counterpart HD-5103B where, as a consequence of 
a recent major merger, an important star formation episode 
took place at redshift ^ 0.3. New stars were born in this 
event, and populated mainly the spheroid, which shows an 
U-shaped colour profile and a mean age profile younger than 
the other galaxies. The distribution of metallicity and a- 
elements for stars in the spheroid shows a clear bimodality. 
Such a bimodality seems to be connected with the differ- 
ent origin of stars in this component. The low- metallicity 
peak could be related to the classical bulge, while the high- 
metallicity peak could be associated with a pseudo-bulge 
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Figure 16. [O/Fe] versus [Fe/H] ratios for F and G stars in the 
thin and thick discs (blue and green dots, respectively) of each 
simulated galaxy. 



tion shows that the thick disc is enhanced in a-elements, but 
not as much as that found for the spheroid. 

Galaxies located in low-density regions of the universe 
have thick discs where, for low-metallicity stars ([Fe/H]< 
—0.4), the abundance of a-elements is clearly enhanced as 
compared to thin disc stars of the same metallicity. When 
fitting the vertical density profiles of RGB stars, we find that 
objects in low-density regions have thick disc scale heights 
Zh that do not chan^^^ignificantly along the disc. According 
to some authors (|lbata et a l. 2009), Zh values roughly con- 
stant with the radial distance favour the scenario in which 
the thickening results from internal (non- merger- driven) in- 
stabilities in the thin disc. In the case of galaxies located 
in high-density regions of the universe, we find instead a 
thick disc scale height that changes with the galactocentric 
distance, as expected in a merger-induced scenario for the 
formation of this component. For these galaxies, the trend of 
[a/Fe] versus [Fe/H] for the thin and thick discs also present 
a different behaviour for [Fe/H]< —0.4. Nevertheless, now 
the number of stars with very low metallicities in the thin 
disc is much smaller than that found in the same component 
of galaxies in low-density regions. Consequently, the trend of 
[a/Fe] versus [Fe/H] for the thin disc of high-density galax- 
ies does not extend to very low metallicities. This fact seems 
also point out a merger-driven scenario as the main mecha- 
nism for the thick disc formation of galaxies in high-density 
regions. Indeed, the lack of thin disc stars with very low 
metallicities (i.e., very old ages) for such galaxies seems to 
favour a formation model in which an early thin disc was 
almost completely disrupted by a violent event giving rise 
to the thick disc. Afterwards, a new stellar thin disc was 
formed from the remaining and already metal-enriched gas. 
In order to discern the role of the various mechanisms of disc 
thickening a more detailed analysis of the galaxy evolution 
at different redshifts evolution is needed. Such an analysis 
will be reported in a future paper. 



formed from instabilities of the inner disc. A detailed anal- 
ysis of the formation mechanism and time evolution of the 
bulge component will be presented in a forthcoming paper 
(Obreja et al. submitted). 

The thin disc appears in our simulations as the youngest 
and most metal-rich component, with roughly solar values 
for metallicity. The median age of stars ranges from 3.8 
to 6.7 Gyr, with a just small fraction older than 10 Gyr. 
This is the only component with a non- negligible star for- 
mation rate 1 M© yr~^) at z = 0. The radial distribu- 
tion of ages and colours show that new stars are forming in 
the thin disc according to the insid e-out paradigm. The U- 
shape of such profiles ind icate fsee jMartme z-Serran o et aD 
I2OO9I : ISanchez -Blazauez et al. 20091) that stars have mi- 
grated across the galaxy reaching high radial distances. 

The thick disc in our simulations contains about 16 per 
cent of the total stellar mass and has intermediate prop- 
erties, in between those of the thin disc and the spheroid. 
From its age and metallicity distributions we find that it is 
populated by older and more metal-poor stars than in the 
thin disc, but younger and more metal-rich stars than in the 
spheroid. The thick disc is not currently forming stars, but 
its history of star formation is quieter than in the spheroid, 
and extends in a long period of time, with a rate comparable 
to that observed in the thin disc. The abundance distribu- 
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